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Abstract 



The cross section for np — > dj is calculated at energies relevant to big- 
bang nucleosynthesis using the recently developed effective field theory that 
describes the two-nucleon sector. The El amplitude is computed up to N 3 LO 
and depends only upon nucleon-nucleon phase shift data. In contrast, the 
Ml contribution is computed up to NLO, and the four-nucleon-one-magnetic- 
photon counterterm that enters is determined by the cross section for cold 
neutron capture. The uncertainty in the calculation for nucleon energies up 
to E ~ 1 MeV is estimated to be < 4%. 
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The radiative capture process np — > d'y is a key reaction in the synthesis of nuclei in 
the early universe. Recently, it has been emphasized by Buries, Nollet, Truran and Turner 
(BNTT) ]]J that the uncertainty in the cross section [[J of np — *> d'y at energies relevant 
for big-bang nucleosynthesis (BBN) is difficult to determine due to the lack of data at low 
energies and the lack of information about theoretical estimates. In determining theoretical 
uncertainties in the abundances of the elements produced in BBN, Smith, Kawano and 
Malaney (SKM) J3] assigned a la error of 5% to the cross section for np —>■ d'y, which was 
also used in the recent analysis of BNTT |L|]. m [0 it was found that this 5% uncertainty 
contributes a significant fraction of the uncertainties in the abundances of elements produced 
in BBN. 

In this work we compute the cross section of np — *> d'y using the recently developed 
techniques of effective field theory in the two-nucleon sector For the energy range 

appropriate for BBN (nucleon energies E N < lMeV) it is appropriate to use the effective 
field theory of only nucleons and photons, as presented in ||, which we denote by EFT(f). 
The cross section of np — > d'y for cold neutrons has been computed in the theory with 
pions |7j and in EFT(/f) up to next-to-leading order (NLO) in the effective field theory 
expansion parameter (s). For cold neutrons the cross section is dominated by Ml-capture 
from the ^So channel via the nucleon isovector magnetic moment. However, in addition to 
the contribution from the effective ranges of both the 1 So and 3 Si channels at NLO, there is a 
contribution from a four-nucleon-one-magnetic-photon interaction with a coefficient, ^L[ M1 \ 
that is not constrained by nucleon-nucleon scattering phase shift data. The observed cross 
section for cold neutrons determines tLi . At higher energies, En ~ 1 MeV, the cross 
section for np — > d'y is dominated by the El-capture of nucleons in a relative P-wave, 3 Pq, 
3 Pi, and 3 P2- I n the energy region relevant to BBN the contributions from both El- and 
Ml-capture are important. 

It is important to emphasize that the results of our calculation look very similar to those 
of effective range theory (ER) [|8],|[| (for a detailed discussion see JT(J). One of the interesting 
results from the recent developments in effective field theory is that ER is seen to reproduce 
the leading orders of any particular amplitude for low energy processes. However, ER fails to 
reproduce the true amplitude at and beyond the order at which there is a contribution from a 
local, multi-nucleon-external- field interaction ||. If multi-nucleon-external- field interactions 
do not enter until very high orders in the EFT(^") expansion, then ER will reproduce the 
observed value to high precision, as is the case for the polarizability of the deuteron fl6l,|TT 



However, if a multi-nucleon-external-field interaction occurs at low orders ER can deviate 
substantially from the true result, as is the case for the the capture of cold neutrons, np — > dy 
|7]], or the deuteron quadrupole moment 0. For this process, the conventional understanding 
of this discrepancy is that important contributions from meson-exchange currents have been 
omitted ||12| . However, in effective field theory, this discrepancy results from the omission of 
four-nucleon-one-magnetic-photon operators that enter at NLO and higher in the expansion. 
For np —>■ d'y at finite but low incident nucleon energy, the two dominant amplitudes, El 
and Ml, behave differently in the effective field theory expansion. In EFT(;f), a four- 
nucleon-one-electric-dipole-photon local operator occurs at N 4 LO, which means that the 
El amplitude can be computed up to N 3 LO with knowledge of only the nucleon-nucleon 
scattering phase shifts. Therefore, this amplitude will look very similar to the expression 
obtained in ER, when a 7/? expansion is performed. In contrast, the Ml amplitude receives 
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a contribution from a four-nucleon-one- magnetic-photon at NLO, and therefore the effective 
field theory result will deviate from that obtained in ER in a significant way. In addition to 
the expressions we obtain for both the El and Ml amplitudes being analytic and compact, 
they are perturbatively close to the true amplitudes for this process, giving a total cross 
section that deviates < 4% over the range of center-of-mass kinetic energies below 1 MeV. 
In our calculation we neglect both isospin violation and relativistic effects, as in the energy 
region of interest both effects are significantly smaller than the uncertainty introduced by 
not computing beyond the order to which we work (relativistic effects are formally NNLO in 
EFT (ft) JJ, but are suppressed by additional factors of m^/Mjy compared to other NNLO 
effects). 

The strong interactions between two nucleons in the 3 Si-channel are determined by the 
lagrange density, up to N 3 LO, 



4 3 *> = -t!™ (N T my (^m) + be? 3 1} [(iv T p 4 iv)t [N T of ] n) + h.c. 

LitfM (N T 0l 2) N)\N T 0l 2) N) - -^U^O^N^^m) + h.c. 
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where Pi is the spin-isospin projector for the 3 Si channel 



Pi = —^(y 2 Oi t 2 



Tr Ps Pi = Ida 



The Galilean invariant derivative operators and are defined by 
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where the covariant derivative is defined to be D = V — ieQA. The superscript on the 
coefficient denotes the number of derivatives in the operator. The expansion parameters of 
EFT(^") are the external momentum involved in the particular process normalized to the 
mass of the pion, Q ~ p/m n . For momenta of order the pion mass or greater this expansion 
will fail to converge. We have not shown the other operator involving six derivatives as it 
does not contribute to np — > <Py at the order to which we are working. In order that the 
deuteron pole is not shifted order-by-order in the Q-expansion, the coefficients appearing in 
eq. ([]]) have expansions in powers in Q , e.g. 



^( 3 Si) 



b 0,-l + °Q,0 t °0,1 ~r 
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The second subscript on each coefficient denotes the powers of Q in the coefficient itself. 
Relating the S-matrix obtained from the lagrange density in eq. ([]]) to that described by the 
effective range expansion, 



|k|cot5 = -7t + \pd(M 2 + l 2 t ) + w 2 (|k| 2 + 7 , 2 ) 2 + 



(5) 
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one can fix most of the coefficients (only one linear combination of €f l} and ^ con- 
tributes to NN scattering) appearing in eq. (0) in terms of 7 t ~ 1 = 4.318946 fm (j t is the 
deuteron binding momentum), pd = 1.764 ± 0.002 fm (the effective range parameter), and 
w 2 = 0.389 fm 3 (the shape parameter) |13| . The neglect of relativistic effects allows us to 
set 7t = 7 = \/MnB where B = 2.224575 MeV is the deuteron binding energy. In addition, 
to the order we are working mixing between the 3 S\ and 3 Di channels does not contribute, 
and so we will not discuss this sector. However, there is a contribution from P-wave final 
state interactions in the El-capture process that enter at N 3 LO. The P-wave interactions 
are described at leading order by the lagrange density 
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where the P-wave operators are 
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(p) 

and Pj is the spin-isospin projector for the isotriplet, spintriplet channel 
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The measured P-wave phase shifts (as given by the Nijmegen phase shift analysis |14[]) fix 
the coefficients appearing in eq. M) to be 
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and 
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Finally, there are interactions with the electromagnetic field that are not simply related 
by gauge invariance to the strong interaction dynamics. The lagrange density describing the 
leading interactions that contribute to np — > d^f is 

£ 2 ,b = [e *L{ M1) (N t P N)\N T P 3 N) B { + h.c. 

+ i e *L[ E 1] (N t 0? a P) N)\N T P a N) E, + h.c. 

- 1 e n E1 \N T 0% P) N)\N T 0?> N) E, + h.c] , (10) 
where E is the electric field, B = V x A is the magnetic field, and where 

P3 = —^2 r 2 r 3 , (11) 



is the projector for the ^So channel. The renormalization group (RG) evolution of *L[ m) 
has been discussed in H, and the evolution of the electric coefficients are determined by 
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These RG equations tell us (by considering the size of the quantities in square brackets at 
/i = run, the matching scale) that the combination tL^ 1 ^ — M^^Oi^i is °f order Q~ l or 
higher, as opposed to the naive counting of Q~ 2 and that is of order Q 

or higher, as opposed to Q~ A . 

The amplitude for np — > g?7 can be written as 



T = eX E1 Ul r 2 r 3 a 2 a ■ e* [d) U v P • + xe X M1 e abc e^ } k b U?t 2 t 3 ct 2 U p 



: (7) 
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where we have not shown amplitudes that contribute much less than 1% to the total capture 
cross section in the energy range of interest, leaving only the isovector El and the isovector 
Ml amplitudes. U n is the neutron two-component spinor and U p is the proton two-component 
spinor. |P| is the magnitude of the momentum of each nucleon in the center of mass frame, 
while k is the photon momentum. The photon polarization vector is e( 7 ), and is the 
deuteron polarization vector. For convenience, we define dimensionless variables X, by 
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In terms of these amplitudes, the total cross section for np — > is 

Ana ( 7 2 + |P| 2 ) 3 
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where we have used nonrelativistic kinematics, as is appropriate for the energy region of 
interest . 

Explicit calculation of |Xei| 2 up to N 3 LO (i.e. Q 1 ) gives 
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The momentum expansion of the theory is made explicit in eq. (|T6T), and it is clear that 
we have captured all terms up to and including (7/m^) 3 , and (\P\/m n ) 3 . Terms that have 



been omitted are of the form (7^^) ~ 0.03, (|P|pd) ~ 0.006 and higher and also relativistic 
corrections of the form (7/M^) 2 ~ 0.002, or (|P|/Mtv) 2 ~ 0.001 and higher, for nucleon 
energies of E ~ 1 MeV. It is interesting to note that contributions from ^ and ifcf 
occur at N 3 LO. However, there are also contributions from the four-nucleon-one-electric- 
photon operators in eq. ([II]) with coefficients that exactly reproduce the renormalization 
scale independent quantities that occur in eq. (|12|). Therefore, these combinations are higher 
order than N 3 LO. The relative contributions from the P-wave final state interactions entering 
at N 3 LO are much smaller than the (7Pd) 3 contribution also entering at N 3 LO. 

The Ml-capture contribution |Xa/i| 2 has been computed up to NLO (i.e. Q 1 ) for the 
capture of cold-neutrons in the theory with pions [0 and EFT(yf) |J. It is straightforward 
to extend these results to the capture of nucleons with non-zero momentum, and at NLO 
we find 
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where Ki is the isovector nucleoli magnetic moment, ai = — 23.714±0.013 fm is the scattering 
length in the ^So channel and tq = 2.73 ± 0.03 fm is the effective range in the ^So channel. 
The constant ^L np is a RG invariant combination of parameters, 
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that must be determined from data. It is simplest to determine ^L np from the cross section 
for cold neutron capture, that is dominated by the Ml matrix element. For incident neutrons 
= 2200 m/s the cross section for capture by protons at rest is measured to 



334.2 ± 0.5 mb |I5[|. The value of tL np required to reproduce a cross section of 



with speed 
be cr expt = 

334.2 mb is ^L np = —4.513 fm 2 (the experimental uncertainty in this measurement introduces 
a negligible uncertainty in our predictions and we have used the isopin averaged value of 
the nucleon mass, = 938.92). The largest uncertainty in the Ml contribution to the 
total cross section is expected to be of the form |P|7p^, which is ~ 10% at an energy of 
1 MeV (corrections of the form ^ 2 p\ have been renormalized away by fitting ^L np to the cold 
neutron capture cross section). Given that for energies above ~ 200 keV the El amplitude 
dominates the cross section, a relatively large uncertainty in the Ml amplitude does not 
lead to a large uncertainty in the total cross section. In fact, we find that by assigning an 
uncertainty of 3% to the El cross section and an uncertainty of the form IPI7P 2 , to the Ml 
contribution, the uncertainty in the total cross section is <• 4% (the maximum uncertainty 
occuring at an energy of ~ 200 keV), where we have added the errors linearly. For energies 
below 8 keV the uncertainty drops to below 1%. 

At LO, the expression in eq. ( pD for the Ml amplitude reproduces the analogous ex- 
pression found in ER |8||| (after the typographical errors in have been corrected ||). 
However, at NLO where the effective range parameters pd and first appear the expression 
differs both qualitatively and quantitatively. The ER amplitude does not correctly describe 
physics at distance scales of order l/m n or shorter. In EFT^), physics at such distance 
scales, beyond the physics of NN scattering alone, is reproduced by the local counterterm 

A comparison between the cross section obtained with effective field theory and the 
numerical values obtained from the on-line nuclear data center || is shown in table |I[ One 
sees that the analytic expressions we have obtained reproduces very well (within a few 
percent) the numerical values of ref. ||. However, at E = 1 keV the EFT(/f) cross section 
is ~ 2.5% lower than the value from 0. The uncertainty in the EFT(yf) value is expected 
to be 7|P|p^ ~ 0.4%, and so it is very unlikely that higher order contributions will bring 
the EFT(tF) value into agreement with the value from ref. |J at this energy. A similar, but 
much weaker statement can be made about the cross section at E — 10 keV. 

A measure of the accuracy of our calculation can be determined by examining the 
deuteron photo-dissociation cross section, which is related to the capture cross section by 

2M N (E J - B) 
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FIG. 1. The cross section for np — > d'j as a function of the center- of-mass kinetic energy E in 
MeV. The dotted curve is the contribution from Ml-capture, the dashed curve is the contribution 
from El capture and the solid curve is the sum of the Ml and El capture cross sections. Both the 
vertical and horizontal axes are logarithmically scaled. 
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Ml (mb) 
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Ml+El (mb) 


(mb) 


1.264 x 10" 8 


334.2 


5.1 x 10" 6 


334.2 W 


332.0 


5.0 x 1(T 4 


1.668 


1.0 x 10~ 3 


1.669(4) 


1.660 


1.0 x 10~ 3 


1.171 


1.42 x 10~ 3 


1.173(4) 


1.193 


5.0 x 10~ 3 


0.496 


3.17 x 10~ 3 


0.499(4) 


0.496 


1.0 x 10~ 2 


0.329 


4.48 x 10~ 3 


0.333(4) 


0.324 


5.0 x 10~ 2 


0.0987 


9.84 x 10~ 3 


0.109(3) 


0.108 


0.100 


0.0501 


0.0136 


0.064(2) 


0.0633 


0.500 


0.00803 


0.0260 


0.034(1) 


0.0345 


1.00 


0.00375 


0.0310 


0.035(1) 


0.0342 



TABLE I. The cross section for np — * dj in millibarns as a function of the nucleon cen- 
ter-of-mass energy, E. The counterterm ^L np is fit to reproduce a cross section of 334.2 mb at an 
incident neutron speed of |v| = 2200 m/s. The fact that this cross section is an input is denoted by 
the asterisk W. The numbers in parenthesis are the uncertainty in the last digit, and are estimated 
by assigning a fractional error of ('ypd) 4 = 0.028 to the El cross section and a fractional error of 
7|P|p^ to the Ml cross section. The last column is the total cross section as extracted from the 
on-line nuclear data center 0. 



where is the incident photon energy and the deuteron is at rest. A comparison between 
the low-energy cross section computed with EFT(/f) and high precision experimental values 
can be seen in table |I[ A detailed and very illuminating discussion of the experiments that 
contributed to these data points can be found in ref . [JTTJ . In addition, a detailed comparison 
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FIG. 2. The cross section for yd — > np near threshold, as a function of the incident photon 
energy in MeV. The solid curve corresponds to the cross section computed in EFT(^f). The two 
dotted curves correspond to the uncertainty in the EFT(/r) calculation as estimated by the method 
described in the text. The two data points with error bars can be found in table |7J. 
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E 1 (MeV) 


Ml (mb) 


El (mb) 


Ml+El (mb) 


(mb) 


2.62 


0.380 


0.866 


1.25(5) 


1.300 ± 0.029 


2.76 


0.327 


1.167 


1.50(6) 


1.474 ± 0.032 


4.45 


0.128 


2.353 


2.48(9) 


2.43 ±0.17 


5.97 


0.084 


2.122 


2.21(8) 


2.162 ±0.099 


6.14 


0.081 


2.084 


2.17(8) 


2.190 ±0.100 


7.25 


0.066 


1.838 


1.90(7) 


1.882 ±0.011 


7.39 


0.065 


1.808 


1.87(7) 


1.840 ±0.150 


7.60 


0.063 


1.764 


1.83(7) 


1.803 ±0.016 


7.64 


0.062 


1.756 


1.82(7) 


1.810 ±0.028 


8.14 


0.058 


1.656 


1.71(6) 


1.800 ±0.130 


8.80 


0.053 


1.534 


1.59(6) 


1.586 ±0.011 


9.00 


0.052 


1.500 


1.55(6) 


1.570 ± 0.036 



TABLE II. The cross section for yd — ► np in millibarns computed in EFT(/f) and the experi- 
mental values taken from pages 78 and 79 of ref. [ 10 1 . The EFT(/r) cross section is comprised of the 
El amplitude computed to order N 3 LO and the Ml amplitude computed to NLO. The theoretical 
uncertainties are estimated in the same way as those in table |[ 



between the predictions of potential models, in particular the Bonn r-space potential, with 
these data can also be found in 
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For the two lowest energy data points the uncertainties 
in the EFT {-ft) calculation are seen to be nearly a factor of two larger than the experimental 
uncertainties at this order. A plot of the break-up cross section in Fig. @, along with 
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0.5 



4 6 8 10 

E y (MeV) 

FIG. 3. The cross section for "yd — > np as a function of the incident photon energy in MeV. The 
solid curve corresponds to the cross section computed in EFT(/f). The two dotted curves correspond 
to the uncertainty in the EFT (ft) calculation as estimated in the text. The data points with error 
bars can be found in table [||. 



the two low-energy data points, clearly shows the need for more data in this low-energy 
region. A few more high precision measurements between ~ 2.5 MeV and ~ 4.0 MeV would 
provide important constraints on the Ml and El amplitudes in the energy region relevant 
to big-bang nucleosynthesis. It does appear that the ~ 3% uncertainty in the El cross 
section that we have estimated to arise from unknown higher order contributions may be 
an over-estimate. Fig. fl3|) hints that a la error of 1% or 2% might be appropriate, but this 
cannot be justified from EFT (ft) alone. Given that the EFT(ft) at N 3 LO reproduces the El 
amplitude very well, and that there are only two precise data points in the energy region 
that is sensitive to the Ml amplitude, it is unlikely that pushing the EFT(ft) computation 
to one higher order would lead to a noticeable difference in cross section. Even at this order 
in the EFT (ft) expansion, we find good agreement with potential model calculations fli"0 |. 

Returning to the np — ► dj capture process, for energies above ~ 300 keV the cross section 
is dominated by El capture and hence the uncertainty in our calculation is essentially the 
uncertainty in the El cross section, ~ 3%. In order to further reduce this uncertainty 
down to ~ 1% a N 4 LO calculation of the El amplitude is required. As we have discussed 
previously, at N 4 LO there is a contribution from a local four-nucleon-one-electric-photon 
interaction that is not constrained by nucleon-nucleon scattering, but could be determined 
by the deuteron photo-disintegration cross section. The appearance of such an operator is 
not restricted to effective field theory as such interactions will also arise in potential model 
calculations of this amplitude, making roughly the same size contribution. An example of 
this is the deuteron quadrupole moment |TE|JT7[| . To describe the lower energy regime to 
higher precision, the Ml amplitude will need to be computed to N 2 LO or higher. At this 
order there will be additional counterterms, beyond ?L np that will need to be determined by 
data. Fig. (|2|) suggests that the existing two data points in this region may not be sufficient 
to achieve this. 

In conclusion, we have examined the radiative capture process np —>■ d*y in the pionless 
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nucleon-nucleon effective field theory. An analytic expression for the cross section in the 
energy region relevant to big-bang nucleosynthesis is presented and is expected to reproduce 
the true cross section at the few percent level. It is simple to relate the rate for np — > 
to the cross section for — > np and our calculation agrees very well with the existing 
data. There is motivation to perform high precision measurements of the deuteron photo- 
disintegration cross section at very low energies, to tightly constrain the Ml contribution to 
the cross section. 



We would like to thank Gautam Rupak and Scott Buries for useful discussions. We 
would also like to thank Baha Balantekin, Wick Haxton and Brad Keister for bringing this 
issue to our attention. This work is supported in part by the U.S. Dept. of Energy under 
Grants No. DE-FG03-97ER41014. 



10 



REFERENCES 



S. Buries, K. M. Nollet, J. W. Truran and M. S. Turner, Phys. Rev. Lett. 82, 4176 
(1999). 

ENDF/B online database at the NNDC Online Data Service, 
[http : / /www . nndc . bnl . gov . 

M. S. Smith, L. H. Kawano and R. A. Malaney, Astrophys. J. Suppl. Ser. 85 219 (1993). 
D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Lett. B 424, 390 (1998); Nucl. Phys. 
B 534, 329 (1998). 

D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Rev. C 59, 617 (1999). 
J. -W. Chen, G. Rupak and M. J. Savage, |tiucl-th/9902056| . 
M. J. Savage, K. A. Scaldeferri, Mark B. Wise, Nucl. Phys. A 652, 273 (1999). 
H. A. Bethe, Phys. Rev. 76, 38 (1949); H. A. Bethe and C. Longmire, Phys. Rev. 77, 
647 (1950). 

H. P. Noyes, Nucl. Phys. 74, 508 (1965). 

Photodisintegration of the Deuteron: A Review of Theory and Experiment, by H. Aren- 
hovel and M. Sanzone, published by Springer- Verlag, 1991, ISBN 3-211-82276-3. 
J.W. Chen, H. W. Griesshammer, M. J. Savage and R. P. Springer, Nucl. Phys. A 644, 
221 (1999);A%d. Phys. A 644, 245 (1999). 
D. O. Riska and G. E. Brown, Phys. Lett. B 38, 193 (1972). 
J.J. de Swart, C.P.F Terheggen and V.G.J. Stoks, |nucl-th/9509032| . 
V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen and J.J. de Swart, Phys. Rev. C49 
(1994) 2950. 

A. E. Cox, S.A.R. Wynchank and C.H. Collie, Nucl. Phys. 74, 497 (1965). 

B. S. Pudliner, V.R. Pandharipande, J. Carlson, S.C. Pieper and R.B. Wiringa, 
[nucl-th/9705009 . 

[17] D. R. Phillips and T. D. Cohen, |nucl-th/9906091 . 



11 



